Schwinger Boson approach to the fully screened Kondo model. 
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We apply the Schwinger boson scheme to the fully screened Kondo model and generalize the 
method to include antiferromagnetic interactions between ions. Our approach captures the Kondo 
crossover from local moment behavior to a Fermi liquid with a non-trivial Wilson ratio. When 
applied to the two impurity model, the mean- field theory describes the "Varma Jones" quantum 
phase transition between a valence bond state and a heavy Fermi liquid. 
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Recently, there has been a broad growth of interest in 
the behavior of heavy electron materials at a magnetic 
quantum critical point 0. Several observations can not 
be understood in terms of the established Moriya Hertz 
theory of quantum phase transitions^, 0, 0K including 
the divergence of the heavy electron masses pj, the near 
linearity of the resistivit yEl B B B and E/T scaling in 
inelastic neutron spectra 10] . The origins of this failure 
are linked to the competition between the Kondo effect 
and antiferromagnetism and may indicate the emergence 
of new kinds of excitation[J [H El El 

In this letter, we show how to unify the Arovas 
Auerbach|lfil] description of quantum spin systems with 
the physics of the Kondo model by using Luttinger Ward 
techniques and a Schwinger boson spin representation. 
Our goal is to develop a large-TV expansion that con- 
tains the physics of antiferromagnetism and the Kondo 
effect in the leading approximation. Traditional pseudo- 
fermion [lil IrH fl8| representations of the spin are ill- 
suited to a description of local moment magnetism. By 
contrast, a Schwinger boson scheme works well for mag- 
netism, but to date, has not been successfully applied to 
the fully screened Kondo model. We use a method co- 
invented by one of us (0P)^3 i n which the Kondo effect 
is captured in a large N Schwinger boson scheme using a 
multichannel Kondo model where the number K = kN 
of screening channels scales extensively with N. By tun- 
ing the number rib of Schwinger bosons from rib < K to 
rib > K, one is able to describe both the overscreened and 
underscreened Kondo models, however, difficulties were 
encountered in the past work that appeared to prevent 
the treatment of the perfectly screened case, rib = K. In 
this letter, we show how these difficulties are overcome. 

Consider the multichannel Kondo impurity model, 

U = E HH va H»« + H i- Hn b ~ 2S), 
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channel index v G spin index a € where 

N = 2j + 1 is even. ipl a = £ k c^ va creates an 
electron in the Wannier state at the origin, where J\f s is 
the number of sites in the lattice. The operator creates 
a Schwinger boson with spin index a € [—j, j] . The local 
spin operator is represented by S a p = b^bp — S a p/N and 
the system is restricted to the physical Hilbert space by 
requiring that rib = Y] a b a b a — IS. The final term in 
H contains a temperature-dependent chemical potential 
A(T) that implements the constraint (rib) — 2S*. We will 
examine the fully screened case 25 = K, taking the large 
TV limit where N —> oo keeping k — K/N fixed. 

We begin by factorizing the interaction in terms of aux- 
illiary Grassman field x„, 



ff ^E^[(^Mxt+H.c.]+^ 
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Following the steps outlined by us in earlier workl2r1 
we now write the Free energy as a Luttinger Ward|2£ 
functional of the one particle Green's functions, 



F[0]=TStr [hi(-g- 



Q- 1 )Q}+ Y[G] (3) 



where Str[A] = TrL4#] — TrL4i?] is the graded (super) 
trace over the Matsubara frequencies, internal quantum 
numbers of the bosonic (B) and fermionic (F) com- 
ponents of A. Go is the bare propagator and Q = 
Diag[G fc , G , G_c], the fully dressed propagator, where 



G b {ivn) 
G x (iuj n ) 



= [iu n + A - E h (iu n )] 1 S crcr >, 

= [G~o(iLj n ) - Y, c (iuj n )\~ 1 8 ac! ,5 ul ,> . (4) 



Here c- creates a conduction electron of momentum k, 



S(iw n ) = Gq 1 — G 1 denotes the corresponding self- 
energies. G c o(iuj n ) — J2k ^ - 1S t ne bare conduction 
electron Green's function. The quantity Y[Q] is the sum 
of all closed-loop two-particle irreducible skeleton Feyn- 
man diagrams. In the large N limit, we take the leading 
0(N) contribution to Y (Fig. [TJ. 

The variation of Y with respect to G generates 
the self-energy 6Y/SQ = S, which yields E c (t) = 
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FIG. 1: Leading contributions to Y[Q\ in 1/N expansion. 
Solid, dashed and wavy lines respectively represent G c , G x 
and Gb ■ Each vertex is associated with a factor i / %/N. Brack- 
eted terms are dropped in the large N limit. 



jjG x (— T)Gfc(r). Since S c is of order 0(1/N) we can 
use the bare conduction propagator G c q inside the self- 
consistent equations for S x (r) = Gf,(r)G c o(— r) and 
Sfe(r) = —kG x (r)G c0 (T). In terms of real frequencies, 
these expressions become 



S x (w — id) 
E b (v - iS) 



^[h F {v)G'^{v)G b {Lu + v) 
h B {v)G'£{v)G* cQ {v - u)] 
k [ — h F {u)[G'/ G {u)G x {v - u) 



+ G'^)G c0 {p-u)] 



(5) 



where G{oj) = G(uj — iS) and the primed variables denote 
the real and imaginary parts of G. We use the notation 
tiB,F(u) — ^±ub,f{^) where ub,f(^) — (e' 3 " =F l) -1 are 
the Bose and Fermi occupation numbers. The constraint 
n b = 25 now becomes / —nz(v)G'l(v) = ^. 

The original work in 0, focussed primarily on the 
case of the overscreened Kondo model, where K > 25. 
The perfectly screened case where K = 25 presented 
two difficulties. First, the phase shift associated with 
the Kondo model is tt/N, which vanishes in the large N 
limit. Second, the requirement that K = 2S appeared ex- 
tremely stringent, the slightest deviation from this condi- 
tion apparantly leading to underscreened or overscreened 
behavior at low temperatures. 

There are two new observations that enable us to now 
avoid these difficulties. First, the perfectly screened case 
where 25 = K is a stable "filled shell" singlet con- 
figuration of the spins. In strong-coupling, this sta- 



ble ground-state corresponds to a singlet "rectangular" 
Young tableau representations of SU(N). [22I [27]]. In 
the gauge theory description of the Kondo model, this 
stability manifests itself as a gap in the Schwinger boson 
and x fermion spectrum. When the chemical potential, 
A lies within this gap, the ground-state Schwinger boson 
occupancy locks into the value n b = K = 25 (Fig. 2a). 

The gap in the gauge particle spectrum has the effect 
of "confining" these excitations, so that the low energy 
physics only involves the elastic scattering of the conduc- 
tion electrons off the Kondo singlet. From this perspec- 
tive, the fully screened Kondo model is a kind of "spinon 
insulator" in this large- N description, while over and un- 
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FIG. 2: (a) Schematic variation of occupancy as a functional 
of boson chemical potential A in the ground-state, (b) Gap in 
spinon spectrum develops in the fully screened Kondo state. 



derscreening develop when the spin-chemical potential is 
in the "valence" and "conduction" bands respectively ( 
Fig. 2 (b)). A second aspect to the problem concerns the 
scattering phase shift. Although the conduction phase 
shift is tt/N, its effect on the thermodynamics is en- 
hanced by the N spin components and the K scattering 
channels, producing an order 0(N x K/N) = 0{N) ef- 
fect in the mean field theory. Moreover, we can identify 
exact Ward identities j2|| which give rise to a sum rule 
relating the scattering phase shift of the conduction elec- 
trons S c = Imln[l — .9 c o(0 — i<5)E c (0 — i<5))] to the phase 
shift S x = Im ln(l + J/<S x (0 — iS)) associated with the 
X fermions, 8 C — -^8 X . The confined nature of the x 
fermions means that S x = n in the ground-state, which 
guarantees that the Friedel sum rule S c — tt/N is satisfied 
in this Schwinger boson scheme. 

In the large N limit the entropy j2|| is given by 



S(T) 
N 



*> f [i mm (-G 6 - 1 ) + G' 6 £' 6 '] + k dnF{w) 



dT 



dT 



Imln(-G^) + G'Z"-G%£' c 



(where the frequency labels oj in the integrand have been suppressed), and 



E c (lj — id) = NY, c (lu — iS) 



dv 



— [h F {v)G'^{v)G h {w + v) + h B {v)G'l{v)G x {v - u)] 
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is the rescaled conduction self-energy. At low temperatures, the gap in the boson and x fermion spectrum means that 
only the conduction electron contribution dominates the entropy, and this is the origin of the Fermi liquid behavior. 
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We can also calculate the local magnetic susceptibility 

XiocCO = -2N J ^h B (u)G' b {u)G' b \u) (8) 

where we have taken the magnetic moment of the local 
impurity to be M = ab^ba, where a = sign(c). Note 
in passing that the dynamic counterpart (S(t)S(O)} van- 
ishes exponentially due to the gap in the bosonic spectral 
functions, the l/t 2 term characteristic of a Fermi liquid 
only appearing at the next order in 1/7V. 




T/Tg 



To see how our method handles magnetic correlations, 
we have applied it to a two-impurity Kondo model, 

k,i> 7 ot 

_ (10) 

where Hk(i) is the Kondo hamiltonian for impurity (i) 
and the antiferromagnetic interaction between the two 
moments is expressed in terms of the boson pair operator 
B12 = J2a &bi a b2-a- H is invariant under spin transfor- 
mations in the symmetry group SP (N) (N-even) 

jT3. We 

now factorize the antiferromagnetic interaction|15j. 

-^Bl 2 B 12 ~,AB 12 +B\ 2 A+^^. (11) 

Boson pairing is associated with the establishment of 
short-range antiferromagnetic correlations. Once A be- 
comes non-zero, the local gauge symmetry is broken, and 
the Schwinger bosons propagate from site to site. In 
this "Higg's phase" the \ fermions also delocalize, giv- 
ing rise to a mobile, charged yet spinless excitation that 
is gapped in the Fermi liquid. Loosely speaking, these 
excitations are mobile Kondo singlets or "holons" . Since 
the paired Schwinger bosons interconvert from particle to 
hole as they move, they only induce holon motion within 
the same sublattice. In the special case of two-impurity 
model, so long as the net coupling between the spins is 
antiferromagnetic, the holons will remain localized. 

Under these assumptions, we can adapt the single im- 
purity equations to the two-impurity model by replacing 



FIG. 3: Showing temperature dependence of (a) impurity 
magnetic susceptibility and (b) specific heat capacity in the 
fully screened Kondo model for k = 0.3, 0.5, 0.7 (with the 
Kondo temperature given by Tk = De~ 2D ^ jK , D being the 
electron bandwidth). Inset calculated Wilson ratio. 

We have numerically solved the self-consistent equa- 
tions JSJ for the self-energy by iteration, imposing the 
constraint at each temperature. Fig. 3. shows the tem- 
perature dependent specific heat coefficient ^ = |y and 
the full magnetic susceptibility x(T)- There is a smooth 
cross-over from local moment behavior at high temper- 
atures x ~ ?1 b(l + n b)/T, to Fermi liquid behavior at 
low temperatures. From a Nozieres-Blandin description 
of the local Fermi liquid |2l) (where channel and charge 
susceptibility vanish), we deduce the Wilson ratio 



W 



xh (i + A) 



X0/70 



1 - 



(9) 



This form is consistent with Bcthe Ansatz results |22j. We 
may also derive this result by applying Luttinger-Ward 
techniques |28j to our model. Our large N approximation 
reproduces the limiting large N behavior of this expres- 
sion, W = 1 + k, in other words, the local Fermi liquid is 
interacting in this particular large N limit. 



Gb{u) — ► G b (uS) = (Gb{u)~ l - \A\ 2 G b (-w)*) 1 (12) 

in the integral equations. We must also impose self- 
consistency A = —Jh{B\2), or 

1 „ f duj , , , 
— = -2 / —h B (v)Im 
Jh J t 



(13) 

We have self-consistently solved the integral equations 
with the modified boson propagator. Using the entropy 
as a guide, we are able to map out the large N phase 
diagram for this model(Fig. 4.). 

We find that the development of A ^ preserves the 
linear temperature dependence of the entropy at low tem- 
peratures, indicating Fermi liquid behavior. However, as 
the Jh increases, the temperature range of Fermi liquid 
behavior collapses towards zero, vanishing at a quantum 
critical point where Jh = J c - F° r Jh > J a Fermi liquid 
behavior re-emerges, but the phase shift 5 X is found to 
have jumped from it to zero, indicating a collapse of the 
Kondo resonance. The entropy develops a finite value at 
the quantum critical point which is numerically identical 
to one half the high temperature entropy of a local mo- 
ment, (iV/2)[(l + 7ib) ln(l + nb) — rib In ?!£,]. Similar behav- 
ior occurs at the "Varma Jones" fixed point [H HJ 
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FIG. 4: Phase diagram for the two-impurity Kondo model 
showing the boundary where boson pairing develops. Color 
coded contours delineate the entropy around the Varma Jones 
fixed point. Black line indicates upper maximum in specific 
heat, blue line, lower maximum in specific heat where cross- 
over into the Fermi liquid takes place. Inset - (a) entropy 
for various values of Tk / Jh (b) showing dependence of 8 X on 
Tk/Jh at a temperature T/T K = 0.02. 

in the N = 2 two-impurity model when the conduction 
band is particle-hole symmetric. We can in fact identify 
two maxima in the specific heat, indicating that as in the 
N = 2 Varma Jones fixed point, the antifcrromagnetic 
coupling generates a second set of screening channels, 
leading to a two-stage quenching process. 

The survival of the Varma Jones fixed point at large N 
in the absence of particle hole symmetry is a consequence 
of the two impurity Friedel sum rule which tells us that 

(S + +S.) = - (14) 

where S± are the even and odd parity scattering phase 
shifts. For N = 2, this condition is satisfied with 5 + = it 
and S- = 0, so it is possible, in the presence of particle- 
hole asymmetry, to cross smoothly from unitary scatter- 
ing off both impurities, to no scattering off either, while 
preserving the sum rule. However, for N > 2, where 
S + + S- < ir, the sum rule can not be satisfied in the 
absence of scattering, so the collapse of the Kondo effect 
must occur via a critical point. 

In conclusion, we have shown that a Schwinger bo- 
son approach to the fully screened Kondo model can be 
naturally extended to incorporate magnetic interactions. 
One of the interesting new elements is the appearance 



of mobile, yet gapped "holon" excitations in the antifer- 
romagnetically correlated Fermi liquid. Future work will 
examine whether these excitations can become gapless at 
a heavy electron quantum critical point, leading to quan- 
tum critical matter with spin-charge decoupling [H Il3]| . 
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